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Abstract We set up a left ring of fractions over a certain ring of boundary problems 
for linear ordinary differential equations. The fraction ring acts naturally on a new 
. module of generalized functions. The latter includes an isomorphic copy of the dif- 

ferential algebra underlying the given ring of boundary problems. Our methodology 
employs noncommutative localization in the theory of integro-differential algebras 
and operators. The resulting structure allows to build a symbolic calculus in the style 
of Heaviside and Mikusinski, but with the added benefit of incorporating boundary 
conditions where the traditional calculi allow only initial conditions. 
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1 Introduction 

General Context. Linear boundary problems are a crucial concern of applied math- 
ematics [15 1.2 50 1 . One might thus expect a rich algebraic theory with symbolic 
algorithms, so as to support the exact solution and manipulation of a suitable class of 
boundary problems. Alas, this is not the case yet. 

There may be two reasons for this. One is that the classical algebraic theory of 
(nonlinear) differential equations, the differential algebra built up by Ritt 11371 and 
Kolchin l24l . does not lend itself easily to boundary conditions: The elements of a 
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differential ring or field are not functions but abstract objects that cannot be evaluated 
at a boundary point. 

The other reason is of a pragmatic nature. If a differential equation can be solved 
at all, one tends to delegate boundary conditions to adhoc postprocessing steps that 
would adapt the integration constants/functions of the "general solution". While this 
may be viable for linear ordinary differential equations (LODEs), the notion of gen- 
eral solution is much less useful in the case of partial differential equations (LPDEs): 
For example, there is not much point in solving the u xx + u yy = per se, but there are 
useful representations of those solutions that satisfy Dirichlet boundary conditions on 
the unit disc (Poisson kernel). So even though we will be dealing only with LODEs 
in this paper, the larger context of LPDEs should be kept in mind. 

In the admittedly modest case of boundary problems for LODEs, an algebraic 
theory was set up in [39|. The decisive step was to expand the well-known structure 
of a differential ring, complementing its derivation by a compatible Baxter opera- 
tor (integral operator). This solves at once two problems: It provides the algebraic 
structure needed for expressing Green's operators (solution operators for boundary 
problems), and it yields an evaluation (a multiplicative functional). The notion of 
evaluation is the key for imposing boundary conditions on the otherwise abstract ob- 
jects of a differential ring; one supplies as many evaluation as needed (in addition to 
the one coming from the Baxter operator). 

In 1 40 1, the basic theory of [39 1 was refined and generalized. Moreover, the au- 
thors have introduced a multiplicative structure on boundary problems that will also 
be crucial for our buildup of the Heaviside calculus. The multiplication of boundary 
problems corresponds to the composition of their Green's operators in reverse order; 
see Equation ©. 

The framework of boundary problems of [39 40 1 was implemented several times: 
The first implementation was coded in MATHEMATICA/TH30REMV as an exter- 
nal package for boundary problems with constant coefficients. This version was su- 
perseded by a new implementation as an internal TH30REMV functor for generic 
integro-differential algebras in [42, 52,(43). Now in the third generation, the most 
recent implementation is the Maple package IntDiffOp [26,27 28 1, which provides 
additional support for singular boundary problems (all previous packages being re- 
stricted to regular boundary problems). 

Heaviside Calculus. Our main concern in this paper is to build a bridge from the 
framework of symbolic boundary problems to the classical Heaviside calculus ETl 
I22II3TI . Since O. Heaviside's main idea, subsequently made rigorous by J. Mikusinski, 
was to treat the differential operator as a symbolic multiplier, it is perhaps natural to 
reflect on possible connections to contemporary symbolic methods for differential 
equations. But why to symbolic boundary problems? 

The answer lies in Heaviside's so-called fundamental formula, the algebraic ana- 
log of a well-known relation for the Laplace transform. If s is the "symbolic mul- 
tiplier" representing the differential operator and / is a suitable function, then one 
has sf = f +/(0) 8q. The point is that s, which is supposed to be invertible in this 
calculus, must somehow "remember" the integration constant that is lost in differ- 
entiation. We will come back to this point in more detail (Section |5). At this point 
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it suffices to say that an appropriate action of s on functions involves an evaluation 
operator. It is used in the Heaviside-Mikusinski calculus for incorporating the initial 
values associated to a differential equation, thus yielding at once the solution of the 
whole initial value problem. It is thus natural — staying in the frame of LODEs — 
to ask if one can build up a more general calculus that would allow to incorporate 
multiple boundary values, given by several evaluation operators. 

The basic ideas of Heaviside and Mikusinski have been vastly extended, specif- 
ically by L. Berg [5] and I. Dimovski fl2l . The latter has also broached a question 
closely related to the one raised above, namely the setup of nonlocal convolutions lfl3l 
l49l and custom-tailored convolutions for boundary problems lfT4l . As far as we know, 
however, there are no direct generalizations of the fundamental formula from one to 
more evaluation operators acting on univariate functions. 

Basic setup. Our own approach is different in many respects. It does not qualify as 
an operational calculus in Dimovski's understanding [ 14], for whom its crucial fea- 
ture is that "operators" and "operands" are merged in a single data structure, whereas 
algebraic analysis keeps the module of operand functions separate from the ring of 
operators. In this sense, we follow the line of algebraic analysis since we distinguish 
the ring of methorious operators (Section |4]i from the module of methorious func- 
tions (Section 01. Moreover, our approach is genuinly algebraic while the hallmark 
of the Heaviside-Mikusinski tradition is an ingenious mix of algebra and analysis. 
Nevertheless, we believe that our setup is close in spirit to the original Heaviside- 
Mikusinski setup, enjoying the following attractive features: 

- Its basic philosophy follows closely in Heaviside's footsteps, making the differ- 
ential operator invertible by "remembering" suitable boundary data. 

- It provides an algebraic structure that accommodates boundary problems for an 
arbitrary number of evaluation points as well as nonlocal conditions. 

- All boundary problems are covered uniformly, so one need not set up custom- 
tailored multiplications (convolutions) for each type of boundary conditions one 
wants to consider. 

- It illuminates the passage from one to several evaluations algebraically: Localiza- 
tion takes place in a noncommutative ring. In contrast, convolution algebras are 
commutative by definition lfl2l . 

- The construction is generic in the coefficient algebra (it works for the class of 
umbral integro-differential algebras). 

From the viewpoint of analysis, however, the class of umbral coefficient algebras is 
rather limited compared to Mikusinski's setup (continuous or even L 2 functions on 
the positive half-axis). While this might be relaxed by suitable limit considerations, 
this is not in the interest of the present paper, where we want to focus on the algebraic 
aspects. 

Our goal in this paper is to build a first bridge between the algebraic theory of 
boundary problems [ 39 , 40'j and Heaviside's tactic of using "symbolic fractions" (lo- 
calization in a suitable convolution ring) for integrating differential equations with 
initial/boundary data. In fact, it gives a new justification (see after PropositionfTOil for 
a notational device initiated in [40|: Since then we have written S3~^ for the Green's 
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operator associated to a boundary problem 38, originally chosen in view of the anti- 
isomorphism (0 mentioned above. 

The work in this paper may also be seen as an answer to an (implicit) question 
originating from |36|. The setting there was restricted to K[x] coefficients and treated 
via Ore algebras. The object of central interest was the integro-differential Weyl alge- 
bra, which has shown to have two important quotient algebras: The ring of localized 
differential operators ^T[d,<? _1 ][j<:] and the usual ring of integro-differential opera- 
tors K[x] [d, f]. The former has d as a two-sided inverse, but no action on K[x]; the 
latter has the action, but J is only a right inverse of d. Now the question arises: Can 
one build up a structure, with an action on K[x] or other coefficient algebras, such that 
the derivation has a two-sided inverse? We show in this paper that the answer is affir- 
mative, provided the derivation is enhanced in a way similar to Heaviside's symbolic 
multiplier s. 

As always in algebra, localizing a ring sheds new light on its structure, especially 
in the noncommutative setting (here the theory of fraction rings is somewhat more 
delicate due to the Ore condition, see Section|4|i. We hope this will also be the case 
for our present construction. But one must bear in mind that the subject has only been 
touched and various issues remain in a preliminary and unsatisfactory state (see the 
Conclusion). Nevertheless, a wealth of new relations has been uncovered, and we are 
confident that they will allow interesting generalizations and refinements. 

Structure of the paper. We start out by summarizing the basic theory of symbolic 
boundary problems, giving special emphasis to their monoid structure (Section 
Our construction is based on a certain kind of boundary conditions, which we have 
called "umbral" because of certain relations to the umbral calculus (Section[3j. After 
this preparation, we tackle the task of localization for an integro-differential algebra 
with an umbral character set (Section|4|i. For the resulting ring of fractions, we con- 
struct a module of functions on which it acts naturally (Section|5]l. We conclude with 
some remarks about possible extensions and generalizations (Section[6]l. 

Notation. All rings are with unit but not necessarily commutative. A domain is a 
(commutative or noncommutative) ring without zero divisors. The zero vector space 
of any dimension will be denoted by O. We use the notation &\ < j^j for indicating 
that a vector space #j is a subspace of a vector space J^j- 



2 The Monoid of Boundary Problems 

Our basic setting is that of fiD . We review the main results here for making our 
present treatment more self-contained and for introducing various pieces of notation 
and terminology in their proper places. We start from the notion of integro-differential 
algebra, a natural generalization of differential algebras that permits the algebraic 
formulation of linear boundary problems. 
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Definition 1 We call d, J) an integro-differential algebra over K if & is a com- 
mutative /f-algebra with /f-linear operations <9 and J such that the three axioms 

(//)'=/, (1) 
{fg)' = f'g+fg', (2) 

(Sf)(Sg')+Kfgy=iSf)g+f(W) (3) 

are satisfied, where . . .' is the usual shorthand notation for d. 

This definition differs from Definition 4 in BP as it uses ([3} as the differential 
Baxter axiom rather than the variant 

Sf8 = fte-Sf'te (4) 

used earlier. In fact, one checks immediately that for commutative & both axioms 
are equivalent (assuming the other axioms). The advantage of (0 is that it is more 
symmetric and that it also sufficient for the noncommutative case, where (|4|i must be 
supplemented by the corresponding dual axiom. Nevertheless, we will from now on 
tacitly assume that all integro-differential algebras are commutative, as in H41\l . 

Incidentally we note the following convention on precedence: The scope of an 
integral sign J covers all factors to its right, unless otherwise specified. Therefore the 
above term jf'jg is to be parsed as J(/' Jg). This helps avoiding the proliferation of 
parentheses in nested integrals. 

The differential Baxter axiom is in general stronger than the following pure Bax- 
ter axiom 

(//)(/*) = ///* + J*//. (5) 

which is the defining axiom for the so-called Rota-Baxter algebras [20,4 44]. (In 
fact, one may relax the conditions on & to include noncommutative algebras over 
arbitrary commutative rings with one. Moreover, one may add a so-called weight 
term for incorporating the discrete setting where J is, for example, the operator of 
partial summation. In that case, one must of course also adapt (0 to account for the 
weight. Remarkably, the differential Baxter axiom © is the same with or without 
weight.) 

Every integro-differential algebra comes with a multiplicative projector, namely 
the evaluation e = 1 — J<9. In fact, the strong Baxter axiom is equivalent to the weak 
one combined with the multiplicativity of e. The presence of the character e leads to 
the direct sum decomposition 

& = (6) 

where = Im(E) = Ker(<9) is the usual ring of constants while J? = Ker(E) = Im( J) 
is called the ideal of initialized functions (one checks immediately that the multiplica- 
tivity of e is equivalent to J 1 being an ideal). So e is the projector onto c & along J* . 
Moreover, both d and J are "^-linear rather than just /^-linear (for d this is of course 
trivial, but for J this is again equivalent to the multiplicativity of e). 

In this paper, we want to restrict ourselves to boundary problems for LODEs. For 
reflecting this property in the (integro) differential structure, in ATI we have called a 
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differential algebra (J^,<9) ordinary if Ker(<9) = K. We will henceforth assume that 
all integro-differential algebras are ordinary in this sense. 

Restricting ourselves to ordinary integro-differential algebras has a number of 
pleasant implications. First of all, the evaluation is now a multiplicative linear oper- 
ator e: J? — > meaning a character (in the sense "multiplicative linear functional 
on an algebra"). Another consequencence, which will be of some importance later, 
is that the polynomials behave as usual. For any integro-differential algebra we 
set x — Jl £ & . Now the pure Baxer axiom ([5]) ensures that x 2 = 2jjl £ J^, and so 
on. It turns out that the elements of K[x] are "really" polynomials, and they satisfy 
the usual differential equations (see before Eq. (6) and also Eq. (8) in ATI for more 
details). 

Proposition 1 Let & be an integro-differential algebra over K, and let K[x] be the 
subalgebra generated by x = Jl. Then K[x] is isomorphic to the univariate polyno- 
mial ring over K. If & is an ordinary integro-differential algebra, we have moreover 



and K[x] < & is an integro-differential subalgebra. 

Example 1 The standard example of an (ordinary) integro-differential algebra is 3? — 
C°°(R) with du = u' in the usual sense and Jm = P u(^ ) d^ . This integro-differential 
algebra contains many important integro-differential subalgebras, for example the an- 
alytic functions C ffl (]R), the exponential polynomials, and of course the polynomial 
ring R[jc]. An important integro-differential subalgebra of C ffl (lR) is formed by the 
holonomic power series [8,45]. 

Given a differential algebra (jF,<9), one may form the ring of differential opera- 
tors J£"[<9]. It is therefore natural to expect a similar ring of integro-differential opera- 
tors &[d, J], which contains differential operators (needed for specifying differential 
equations) as well as integral operators (needed for specifying Green's operators). 
The evaluation e can be used for writing a condition like 2m' (0) — 3(0) = in the 
form (2e<3 — 3e) u = 0. But for boundary problems one usually needs more than one 
evaluation (see Example[2]below). 

For algebraizing general boundary conditions, we start from a character set 0, 
meaning each <p 6 is a multiplicative linear functional just as e is. In fact, we will 
always assume E£$. Based on the characters of 0, we can from arbitrary Stieltjes 
conditions [41, Def. 14]; their normal form is 



with a^j £ K and /«, £ & almost all zero. The summands with J make up the so- 
called global part, those without the local part. The order of a boundary condition 
is the largest i such that a, ^ in the normal form {§), In the standard setting of 
ExampleQ] a typical Stieltjes condition like 



Ker(a") = [M,...,j t »- 1 ], 



(7) 




(8) 
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is encoded as j8(«) =0 with j3 = Eo<9 2 — 3e_i + 7ei +ej Jx 2 — Ei Je^ + E-i Je A , where 
we have written E a for the standard character / n> f(a). Henceforth we shall always 
consider Example Q] with the character set = {e q | a G R}. 

For each character set one may then build up the corresponding operator 
ring ^<p[d, J]. The details are given in [41, §3]; at this point it suffices to note that 
the operator ring is given by the quotient of a free algebra modulo a certain ideal of 
relations like the Leibniz rule, and also Q as well as (0. The resulting operator ring 
decomposes nicely into a direct sum; see HT1 Prop. 17] for a proof. 

Proposition 2 For any integro-differential algebra (J?",<3, J) and character set 0, 
we have the decomposition ^<p[d , J] = jF[<9] + J^"[J] + (0). 

Here (0) is the two-sided ideal generated by <£> in the operator ring /? = j?,j>[<9, J"]. 
We refer to its elements as boundary operators since they may also be described in 
terms of Stieltjes conditions: One checks first that the collection of Stieltjes condi- 
tions can be characterized as the right ideal |<J>) = 0R. It turns out that (0) is then the 
left ^-module generated by \0), so every boundary operator can be written as /ij8] + 

h/n/3„for some functions f\ /„ E & and Stieljes conditions J3i, ... ,j3„ G (4>). 

In particular, every Stieltjes condition is also a boundary operator. 

The decomposition J^#[d, J] = J^[<3] + ^[J] + (<£) reflects the basic needs for 
an algebraic formulation of boundary problems: We need J^[<3] for the differential 
equation, (0) for the boundary conditions, and J^[J] for the solution (Green's oper- 
ators). Let us see how these ingredients look like for the simplest possible boundary 
problem in the standard setting of Example Q] 

Example 2 The boundary problem 



k(0)=h(1) = 



can be encoded in the standard setting J? = C°° (R) by the differential operator T = 
d 2 G J^[<9] and the boundary conditions Eo = e, Ei G (0). Here we can choose = 
{eo,Ei} or any character set containing that. The Green's operator can be written 
as G = Ax + xB — xAx — xBx G &[[], where we use the standard abbreviations A = J 
and — B = (1 — e i ) J. Note that A is the integral from to x, and B the integral from x 
to 1. In the L 2 setting, A and B are adjoint operators. 

The algebraic description in terms of the two functionals Eo and Ei contains some 
arbitrariness since we can clearly form any linear combination of the two functionals 
without changing the solution operator G. In general, a boundary space is a finite- 
dimensional subspace of J^"* generated by Stieltjes conditions. So the boundary space 
of Example|2]is [eo, Ei]. The lattice of all boundary spaces will be denoted by K0. 

We will usually restrict the coefficient functions of the differential operators to a 
differential subalgebra <S < & so as to ensure solutions in the ambient algebra 
This means we require that S be saturated for & in the sense of ATI Def. 18]. 
Whenever an integro-differential algebra (^,d, j) is specified, it is understood that 
a saturated coefficient algebra S is set aside. 
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Definition 2 A boundary problem is a pair (T,8$) consisting of a monic differential 
operator T £ $[9] and a boundary space 8$ G K<P. 

So the boundary problem in Example[2]is given by (<9 2 , [eo, Ei]). Conversely, we 
can think of a boundary problem (T, 8$) as finding a solution «g<? such that 



Tu=f 



for an arbitrary forcing function / g J^. Of course it suffices that j3 ranges over any K- 
basis of ^£f; this is how the boundary conditions are normally given in the first place. 
The traditional formulation in terms of u and / is more intuitive, but it conceals the 
fact that we are really working in the operator ring J£"$ [<?,/]■ 

Clearly every differential operator T has a certain order written as ord T, which 
as in Example|2]usually coincides with the number of given boundary conditions (the 
dimension of 83). In fact, this a necessary — but in general not sufficient — condition 
for the boundary problem to be regular in the sense that there is a unique solution u G 
J£" for every given forcing function / G In turns out that this is equivalent to the 
following formulation in terms of the relevant spaces. 

Definition 3 A boundary problem (T,8§) is called regular if Ker(r) + 8§ L = 
and singular otherwise. 

Here we have written for the orthogonal, meaning the space of all / G & 
such that [i (/) =0 for all /3 G 38. In other words, 8§ is the space of admissible 
functions: those that satisfy the given boundary conditions. The orthogonal, together 
with the analogous notion Sf - 1 for subspaces £f < establishes a Galois connection 
between and its dual . This is similar to the situation in algebraic geometry 
with its Galois connection between affine varieties and radical ideals. See [41 , §5] for 
more details. 

The above definition of regularity is not suitable for algorithmic purposes. How- 
ever, it turns out to be equivalent to the following explicit regularity test [35, Prop. 6.1], 
which is well-known in the special case of local boundary conditions ll23l p. 184]. 

Lemma 1 Let (T,83) be a boundary problem over & and choose bases f5\ , . . . , /3« 
for 83 and u\,...,u n for Ker(T). Then (T,3S) is regular iff ord T — A\m8§ and the 
evaluation matrix J3(m) = [j3;(z<y)] G K" x " is regular. 

Regular boundary problems are exactly those that have a Green 's operator in the 
classical sense (see ll27l on generalized Green's operators for singular boundary prob- 
lems). In our algebraic setting, Green's operators can be defined within the operator 
ring &&[d, J]. 

Definition 4 Let (T,8§) be a regular boundary problem over an integro-differential 
algebra Then G G J^<p [d , J] is called its Green 's operator if TG = 1 and Im(G) = 
8S L . 
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The two conditions for Green's operators are just that u = Gf satisfies the dif- 
ferential equation Tu — f as well as all boundary conditions j3 (u) = for /3 £ 33. 
Since we have assumed a coefficient algebra § saturated for one can prove ATI 
Thm. 26] that every regular boundary problem has a Green's operator in ^#[5, J]. 
In a leap of faith, we write G = (T, 3S)~ l just as in ATI p. 533]. The relation to actual 
inverses will become clear in PropositionfTOl 

Apart from the localization, there is another good reason for the notation G = 
(T, 33) ~ 1 . It turns out that boundary problems can be multiplied in such a way that 

((T u 33 l )(T 2 ,33 2 )y 1 =(T 2) 33 2 )- 1 (T 1 ,3$ l )- 1 (9) 

is satisfied. For that purpose, we have defined the multiplication of boundary prob- 
lems in the fashion of a semi-direct product by 

(71,^0(72,^2) - {T 1 T 2 ,M 1 T 2 + M 2 ). (10) 

It is easy to check that one obtains a monoid S[d] x K<P in this way, with (1,0) 
as the neutral element. Furthermore, it turns out that the regular boundary problems 
form a submonoid £[d}® C £\d\ X K<P, and © means that {T,33) M> (J,3§Y X is an 
anti-isomorphism from $ to the multiplicative monoid of Green's operators. For 
a proof of these facts, we refer to BT1 Prop. 27]. 

The monoid algebra K£[d\& will be called the ring of boundary problems (we 
suppress the qualification "regular" since we will not investigate the singular case in 
the frame of this paper). Its elements will supply the numerators in the localization to 
be constructed below (Section|4|i. 

At this point we should also note that a regular boundary problem (T,33) need 
not be well-posed in the sense that its Green's operator G is continuous in the stan- 
dard setting of ExampleQ] For example, consider the regular boundary problem (d — 
l,E d 2 )or 

u' — u = f 
u"(0) = 

in traditional notation. One checks immediately that G = e*\e~ x — e*Eo — <? x Eo<3 E 
[d, [] is its Green's operator, meaning the general solution is given by 

u(x) = r o e*Sf(t;)dt;-(f(0)+f'(0))eX. 

Clearly, this Green's operator G is not continuous, at least not if one endows C°°(]R) 
with the usual topology induced by C(R). Following Hadamard, a well-posed bound- 
ary problem lfl6l p. 86] must be regular (meaning the solution u exists and is unique 
for each given /) as well as stable (meaning the solution u depends continuously on 
/). In the example above, the source of the instability is clear — we have imposed a 
second-order condition for a first-order differential equation. Going back to the al- 
gebraic setting then, we call a boundary problem (T,3§) well-posed if it is regular 
and 33 can be generated by Stieltjes conditions whose order is smaller than that of T; 
otherwise (T,33) is called ill-posed. As one easily checks, the well-posed boundary 
problems form a submonoid of <t>\d]$. 
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The Green's operator G of a regular boundary problem (T, 3$) of order n factors 
naturally as 

G=(l-P)T°, (11) 

where T° is the so-called fundamental right inverse of T, defined as the Green's 
operator of the initial value problem (T, [e,e<3, . . . , Ed" -1 ], and P is the projector 
onto Ker(r) along 3$ L . The fundamental right inverse T exists since initial value 
problems are always solvable [41 , Prop. 22]. We may thus think of (fTTb as a two step 
process: Using T° one solves first the initial value problem (which is usually much 
easier), then one incorporates the boundary conditions via the projector P. 

The submonoid relation C <?[<?] x K& means that multiplying regular 

problems leads to regular problems again. It is interesting to note that the converse is 
also true, provided the order condition of Lemma [TJ holds. 

Lemma 2 Let (Ti,3§\), {T2t8$2) be boundary problems over an integro-differential 
algebra J? with ord7\ = A\m.3§\ and ord^ = d\m.3§2- Then (Ti,3§\) and (T2, 8$%) 
are regular whenever (Tj , 38\ ) • (J2, 3§2) & 

Proof Let us write m = ord T\ = &\m.3§\ and n = ordT? = dim 3§2- Choose fundamen- 
tal systems . . ,/ m 6 & for T\ andgi, . . . ,g n G & for T^- Take ^-bases j8i,... ,j8 m 
of 9»\ < and y u ■ ■ ■ , 7« of ^2 < Then r 2 °/i , . . . , r^/m.Si , • • • ,g n is a fun- 
damental system for T\ T2, and a /^-basis of 3§\ T2 + 3§2 is given by P1T2,..., fimTi, 
Yi, • • ■ , Yn- The latter fact follows since the sum 3§\ T2 + 3§2 in <TT~0b is always direct IT351 
Prop. 3.2]. By Lemma[TJfhe regularity of the boundary problem (TiT^^iT^ + 3§2) 
means that its evaluation matrix 

(m)(i?f) mM\ ( Pi/) \ 

I Y(Tff) Y(g) )-\Y{T?f)Y{g)) 

is regular. But this is only possible if both diagonal blocks j3 (/) and y(g) are regular, 
and these are just the evaluation matrices of (T\,3§\) and (T2,3§2)- 

For practical applications, one is not so much interested in multiplying bound- 
ary problems — thus increasing their order — as to factor them into lower-order prob- 
lems. Interestingly, factorization will also be instrumental in the localization process 
(Lemma [Toi l. We repeat here as the main result the so-called Factorization Theo- 
rem ETJ Thm. 32]. 

Theorem 1 Given a regular boundary problem (T,3§) G <^[d]#, every factoriza- 
tion T = T1T2 of the differential operator lifts to a factorization (T,3§) = {T\,3§\) ■ 
{T 2 ,3§ 2 ) with (J\,88\),{T2,3§2) regular and 3§ 2 < @. 

We note that the right factor may be chosen to be an arbitrary subspace of 38 as 
long as regularity holds for (T2, £%2)\ in particular one can always choose initial con- 
ditions for £82 t° ensure regularity. In contrast, 38 \ is determined by the choice of T\ 
and T2 alone. This is the content of the following uniqueness result [41 Prop. 31], 
called the Division Lemma for boundary problems (see below for the notion of sub- 
problem). 
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Lemma 3 Given a regular boundary problem (T,3S) € <?[d]tf> and any factorization 
T = T1T2 of the differential operator, there is a unique boundary problem (T\,£3\) 
such that for any regular subproblem (J2^i) < (T, 38) we have the lifted factoriza- 
tion (T\,S§\) ■ {T2,S§2) — {T,33). 

The notion of subproblem was defined (under the name "right factor", which is 
avoided here due to the ambiguity in cases like Example IT) for regular boundary 
problems in BTl Def. 29], but it can be extended naturally to all boundary problems. 
Hence a subproblem of a boundary problem (T,3§) £ k K<P is any (72,^2) S 
S[d] k K<P with T2 a right divisor of T and S$2 a subspace of 38. We denote this 
relation by (T 2 ,& 2 ) < (T,@). 

We conclude this section with a few remarks on algorithmic issues. While at the 
present stage of our work we do not aim at a computational realization of the local- 
ization, we see this nevertheless as a mid-term goal. In fact, most of the results in this 
paper are algorithmic, with two exceptions: (1) The Ore condition is needed in (o[d] 
for monic operators; this is not supported by customary packages but appears to be 
adaptable (see the remarks after Proposition|6]l. (2) One would need an algorithm for 
determining the monomials required in Definition [5] Once these two gaps are filled, 
one should be able to compute in the localization since all the constructions reviewed 
in this section are themselves algorithmic [41 1: 

- The operator ring ^0 [d , J] is defined as a quotient by a certain ideal of relations. 
This ideal is described by a noncommutative Grobner basis for an ideal with 
infinitely many generators (on another view: a noetherian and confluent rewrite 
system). 

- If the ground algebra & and the character set are computable, the normal 
forms of ,^&[d, J] are computable as well. Typically, & is an algorithmic frag- 
ment of C°°(IR) like the exponential polynomials l39l p. 176], while <£ consists 
of finitely many point evaluations / M> f(a). 

- The Green 's operator of a regular boundary problem (and hence the solutions for 
arbitrary forcing functions up to quadratures) can be computed as long as one has 
a fundamental system for the underlying homogenous differential equation. For 
the latter, one can in principle rely on the vast body of results from differential 
Galois theory 041 . 

- The multiplication and factorization of boundary problems can be computed as 
long as the same is true for the constituent differential operators. In this way 
one can often solve higher-order boundary problems by decomposing them into 
smaller factors such as in a recent application to actuarial mathematics |3|. For 
the algorithmic theory of factoring linear differential operators, we refer to ifTHl 

HSE3I. 

As detailed in the Introduction, the algorithms addressed above are implemented in 
Mathematica and Maple packages. 

3 Umbral Boundary Conditions 

We have seen that for a given integro-differential algebra & and character set 0, 
the natural choice of boundary conditions is the Stieltjes conditions, defined in Pffl 
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as the right ideal \<P) of J^[<3, J] and characterized by their normal forms (0. The 
choice of \<P) is vindicated by TheoremQ] which asserts that Stieltjes conditions are 
sufficient for describing arbitrary factor problems (and they are clearly also sufficient 
for multiplying boundary problems). In fact, Stieltjes conditions with nonzero global 
part appear in the left factor problem even if the original boundary problem has only 
local conditions; see e.g. Example 33 in [41 1. 

For the purpose of localization we must isolate a subclass of Stieltjes conditions. 
First of all, we must distinguish carefully between the zero condition £ \<P) and 
degenerate boundary conditions that only act as zero: We call a Stieltjes condition j3 
degenerate if <p(f) = for all / £ In the standard setting of ExampleQ] there are 
plenty of degenerate boundary conditions. If / is a bump function supported, say, on 
an interval disjoint from [0, 1], the Stieltjes condition Ei jf — j^f is clearly degener- 
ate. In contrast, the integro-differential subalgebra of analytic functions C m (M) does 
not have degenerate global conditions of this form as we shall see below (Example|4|. 

The subclass of Stieltjes conditions to be chosen must be so as to ensure a suffi- 
cient supply of regular boundary problems. In particular, we shall find ourselves in the 
situation of embedding a singular boundary problem (T, SS) into a surrounding regu- 
lar problem (LemmafTOb. Hence we must enlarge T E <?[<9] by a suitable differential 
operator T € S[d] and the boundary space 38 by new boundary conditions j3i ,j32, . . . 
in such a manner that the resulting evaluation matrix is regular (Lemma[T]i. 

The key for solving this problem is found in the fortunate fact every integro- 
differential algebra & contains the polynomial ring (Proposition [TJ, so the initial 
value problem 

(,?",[£,£,?,..., Ed"- 1 ]) 

is a natural choice for extending the given boundary problem (T, 33). Then the mono- 
mials 1 , xj 1 ! , . . . , x"~ 1 / (n — 1 ) ! are a fundamental system by Proposition Q] and the 
corresponding evaluation matrix is /„. But this is not enough for embedding singular 
problems since the evaluation matrix of the surrounding problem involves combin- 
ing the given boundary conditions j3 € 53 with the monomials of the fundamental 
system (0. Indeed, the crucial step in the proof of Lemma [TOl needs that j3(x m ) ^ 
for some monomial x m G K[x\. This is the motivation for the following definition. 

Deflnition5 A Stieltjes condition j5 £ ^$[<3, J] is called umbral if j3(x m ) ^ for 
some monomial x" 1 6 K[x]. Furthermore, we call <P an umbral character set if every 
nondegenerate Stieltjes condition is umbral. 

The reason for our terminology is that there is an interesting link between Stiel- 
jes conditions j3 € -i^$[<3, J] with j3(x m ) ^ and the umbral calculus. Indeed, every 
umbral condition defines a nontrivial shift-invariant operator (Proposition^. This is 
clear for the local parts q>d" but needs some justification for global terms <p J/ with 
arbitrary / € ^ D K[x]. The point is that one may apply the following special case 
of the well-known "antiderivative Leibniz rule". (In the general setting of analytic 
functions one writes J/g as an infinite series of iterated integrals of / times iter- 
ated derivatives of g. Here we have g — x" so that only finitely many derivatives are 
nonzero and the series terminates.) 
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Lemma 4 In any integro-differential algebra (^,d,J), we have the formula 

J/x"=f(-l) fc nV-V ( - fc - 1) (12) 

£=0 

for all f G Here f~ k (k > 0) is defined by = / and f^ k ~ l) = J f { ~ k) . 

Proof This is essentially the special case m = 1 of Equation (14) in [36], but the proof 
is short enough to present here. We use induction over n. The base case n = is clear, 
so assume (TT2b for n > 0; we prove it for n + l. Since x" +1 = (n+l) JV" we have 

J/x" +1 = (n + 1) J/Jx" = - (n + 1) J/^'V 

by the pure Baxter axiom ((5). Using the induction hypothesis, the right summand 
becomes 

"f^(-l) k {n + l)t x n+l - k f(- k -V 

k=l 

after an index transformation. Adding the extra term p l >x n+l means extending the 
summation to k = 0, . . . ,n + 1, which yields (TT~2T > for n + 1 . 

Lemma 5 Lef j3 = <pj/ £>e a global condition in J^$[<9, J]. Then we have /3 = <j?/3 
as a functional K[x] — > K, where 

P = f,b k d k : K[x]^K[x] (13) 

k=0 

is a shift-invariant operator with expansion coefficients bk — ( — l) k (p(f(~ k ~^). 

Proof By a well-known result of the umbral calculus [11, Thm. 2.1.7], incidentally 
also presented in [38, Prop. 92], the operator j3 is shift-invariant. For seeing that j3 = 
<pj3, we apply Lemma|4]to obtain j3(x") = Y*kbk<?(n-x n ~ k ) = Y.kbk<pd k (x"). 

The operator expansion (fT3l allows us to associate a shift invariant operator with 
a Stieltjes condition of the special form j3 = <pjf. But this association generalizes 
immediately to arbitrary Stieltjes conditions since the local conditions are unprob- 
lematic. 

Proposition 3 Let j3 be a Stieltjes condition in J^"$ [d , J]. Then there is an associated 
shift-invariant operator 

j3 = £M*: K[x]^K[x] (14) 

k=0 

with coefficients = j8 (x k /k\) such that /3 = E j3. Clearly, the associated operator j3 
is nonzero iff /3 is an umbral condition. 
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Proof Let Sq, : K[x] ->K[x] be the shift operator fix) M>/(x + <p) with (p = (p(x)£K. 
Using the normal form ([H) we can write the given boundary condition as 

<pe<P <pe<P 

where Ty = Y.i a <p.i^' is clearly shift-invariant while jS^ is the shift-invariant operator 
corresponding to /3<p = (p$f<p according to Lemma [5] Then the terms SyiTy + fi^) 
in the second sum are shift-invariant, hence we have j3 = Ej3 with the shift-invariant 
operator 

and the formula bk = P(x^/kl) = Ef}(x k /kl) for its expansion coefficients follows 
from the general result referred to in the proof of Lemma|5] 

Umbral boundary conditions appear to be abundant, at least in the cases most cru- 
cial to us, especially in the C°° setting (containing many important integro-differential 
subalgebras — notably the analytic and holonomic functions as well as the exponential 
polynomials). 

Proposition 4 In the standard setting of Example\l\ the point evaluations form an 
umbral character set. 

Proof Consider an arbitrary nondegenerate Stieltjes condition 

<pe<£i'=0 <pe<t> 

so that there is a function with /3 (u) / 0. Let be the maximum of | <p | for all <p 

with a,pj ^ or fm ^ 0. We consider now the Banach space C (K) on the compact 
interval K= [—R,R]. Its norm ||-||*is given by 

11/11* = I> (<) II- 
i=0 

for all feC k (K). There is a little known generalization of the Weierstrass approxima- 
tion theorem due to Nachbin [32], which asserts (in a simple special case) that R[jc] 
is dense in C (K). Hence we may choose a polynomial sequence p„ that converges 
to « e J < C k iK) in the C k topology. 

One checks immediately that j3 : C k iK) — > R is a continous functional with re- 
spect to the C norm. In detail, one has /3 («) | < C||«||^ with operator-norm bound 

c= EL^+LIpIIIaii-. 

<pG<£l'=0 <pG<J> 

Therefore we have /3 (/?„) — ► j3 (m) 7^ 0. This is impossible if j3 vanishes on all of R[jc]. 
Hence there is some p G M[x] with /3(j;>) 7^ 0. Clearly there exists a smallest mono- 
mial x™ in p with /3 (V") 7^ 0. 
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Not every character set is umbral, though. It is natural to introduce the following 
necessary conditions for a character set <£> to be umbral (we write <p for the canonical 
value (p(x) £ K of a character <p S 0): 

1 . The character set <J> must clearly be separative in the sense that (f> = X implies (p = 
X- Otherwise j3 = <p — J is a nonzero Stieltjes condition with j3 (^[jc]) = 0. 

2. Every character must be complete in the sense that every global condition 
of the form j3 = <p J/ is umbral whenever it is nondegenerate. This may also be 
expressed as / K[x] => f = 0. Here orthogonality refers to the nondegenerate 
bilinear form (f\g) = <?//<?■ If this bilinear form is positive definite with K = H 
or K = (D, we have a pre-Hilbert space (j^, (-|-)m). Following the terminology 
of V.24], completeness of <p then means that any <p-orthonormal basis of K[x] 
is complete in (-|-),p). 

These conditions are most likely not sufficient. At the moment we cannot give a coun- 
terexample for corroborating this claim. But it is intuitively clear that the complete- 
ness properties associated with each <p separately cannot prevent linear dependencies 
between the actions of distinct global conditions <p J/. It remains an interesting task 
to formulate stronger conditions on <P that ensure an umbral character set in a natural 
way. 

Example 3 As an example of a nonseparative character set, consider the exponen- 
tial polynomials ^T[x,e x ] with the nonstandard character <p defined by (p(e x ) = 1 
and <p(x") = 1 for all n, effectively mixing evaluation at 1 € K for monomials with 
evaluation at 6 K for the exponential. Choosing K = R for simplicity, let e be the 
honest evaluation / h-»-/(1). Then clearly <p and e coincide on K[x] while they are in 
fact distinct characters on A^e*] since (p(e x ) = l^e = E(e x ). 

This shows that separativity is necessary but not sufficient. Nevertheless, it en- 
sures that all local boundary conditions are indeed umbral. 

Proposition 5 Let <P be a separative character set for an integro-differential alge- 
bra d, f). Then every local boundary condition is umbral. 

Proof Using (O, every local boundary condition can be written in the form 

P = Ha ij(Pi d'- 1 . 

1=1 j=\ 

Assuming f}(K[x]) = 0, we show that all coefficients a, 7 are zero. Collecting them in 
the vector 

a = (01,1,01,2, • • -,a\.s', . . . ; a r \,a r 2, ■ ■ -,a v ) T 6 K n , 

we obtain n = rs linear homogeneous equations for the unkowns ay by applying j3 
to the monomials l,x 7 x 2 /2\, . . . ,x"~ l j(n — 1)!. The matrix of the corresponding sys- 
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tem Ma = is given by M = (M„ s (<p\ ) ■ ■ -M m ((p r )) G A" x " with blocks 
/ 1 \ 



M„ s {x) 



x 

2 



1 



1 



X S X 



— 



(*-!)! 



\(n-l)! (n-2)\ 



X 

£ 

2 



1 

X 



Jt"- 



(n-i+1)! (n-s) 



n-s)! / 



G AW 



(15) 



Since <J> is separative, the values <pi , . . . , <p r are mutually distinct. Then we may apply 
the subsequent Lemma [6] to obtain detM ^ 0, hence Ma = has the unique solu- 
tion a = 0. 



Lemma 6 77ie determinant of M(x) = (M ns (x\) : . . . ,M ns (x r )) G A[xj 
« = ant/ blocks H5i , is given by 

detM(x) = V (rf 2 sf(s - l) r /sf(« - 1), 



with 



(16) 



w/iere V(r) = j | (xj — Xi) is the r x r Vandennonde determinant and sf(/) = 

l<i<j<r 

1!2! • • • i\ denotes the superfactorial. 

Proof The determinant is a special case of |17 Thm. 1.1] or |29, Thm. 20], apart 
from the linear factor 



in k-th row of M(x). Thus we have 



s-l 



<=U=i 



detM(x) = V(rf x ( nrii ! ) ( 11 ^7 ) - W sf(,- l)7sf(»- 1) 



*=o 



which is ( fTBT l. 



Let us now turn to completeness. In the analysis setting, every <p is complete 
since we know already the stronger result that every Stieltjes condition is umbral 
(Proposition|4]i. It is nevertheless instructive to have a closer look at this case. 

Example 4 Consider the point evaluation <p = E a in the standard setting of Exam- 
ple Q] Here we get an inner product 

which can be extended to all continuous functions f,g on [0,a] and so gives rise to 
the pre-Hilbert Hausdorff space (C[0,a], (-|-)<p)- It is well-known that (x m ) m en<s is a 
complete sequence in this space JT] V.24]. We may assume a — 1 by a scale transfor- 
mation. By the usual Gram-Schmidt process, the monomials x m can be transformed 
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to an orthonormal basis (e n ) of C[0, 1], which consists in this case of the Legendre 
polynomials 

*=4^f (**-ir. 

2"n\ dx" y ' 

In every pre-Hilbert Hausdorff space E, an orthonormal basis (e m ) has the following 
well-known property |7, Prop. V.2.5]: If / £ E is nonzero, we have (e m \f) ^ for 
some m £ N. Applying this to E = C[0, 1] D C°°[0, 1] and fixing any / £ C°°[0, 1] C 
C°°(R), we see that <p is complete in the sense defined above. For if (^ n \f)q> = for 
all m £ N, then / vanishes on [0, 1] and hence <pj/ is degenerate. 

The above argument shows in fact that every global condition <pj/ over C°°(R) 
is umbral or degenerate. Working over C ffl (R), we can exclude the degenerate case. 
Indeed, the identity theorem of complex analysis ensures that / = whenever / van- 
ishes on [0,1]. This is in stark contrast to the C°° case, which has plenty of degenerate 
global conditions <p J"/ as we have observed at the beginning of Section[3] 

It becomes clear from the above example that completeness is really an analytic 
property of some sort. It is therefore not surprising that one can construct "purely 
algebraic" examples lacking this property. 

Example 5 To give an example of an incomplete character on & = K[x, e x ] with K = 
R, we proceed similar to Example [3] Define cp(x") — for n > and (p(e x ) = e, 
<p(l) = 1. Then clearly /3 = <j?Jl is nondegenerate since f5(e x ) = (pje x = (p(e x — 1) = 
e - 1 ^ 0. But we have j3 (x m ) = <p JV" = <p(^[) = for all m £ N. 

As noted above, separativity and completeness are most likely not strong enough 
to ensure an umbral character set. Note that umbrality of implies that two Stieltjes 
conditions in J^*[<2, J] are linearly independent on K[x] whenever they are linearly 
independent on J?. If one of the two conditions is local and the other global, it is 
reasonable to expect this property to follow from completeness. This expections is 
fulfilled. 

Lemma 7 Let (J^,<3, J) be an integro-differential algebra with a complete charac- 
ter (p. Then a nondegenerate global condition (pjf never coincides on K[x] with any 
local condition based on (p. 

Proof Assume j3 = (pjf coincides on K[x] with ao<p + a\(pd H ^a s (pd s , where s 

is chosen minimal. Then the umbral expansion (Tol l breaks off at k = s, and we 
have <p(/ (- * -1) ) = for k>s. For / = f- s - l \ we get (p(/ ( "^ 1) ) = for all k £ N, 
so the condition <pj / is degenerate since <p is complete. For showing that this can- 
not happen, it suffices to prove that all <pJ/'~ A ' are nondegenerate whenever q>jf 
is. We use induction on k £ N. The base case k = is trivial, so assume the claim 
for fixed k. By the induction hypothesis we can choose g £ J£" with <p[p-~®g = 1. 
Using (|4| with g in place of / and f(~^ in place of g, we have 

!f^ g =f^ g -Sf^ g ', 
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so our choice of g entails 1 = <p(f'~ k ~^ g) — (pjf^ kl 'g'- If the first summand on 
the right-hand side is different from 1, then g' witnesses to <pjf(~ k ~ 1 ^ being nonde- 
generate, and the induction is complete. Otherwise <p{f^ k ^ 1 ^) must be nonzero, and 
we may use the special case (J7z) 2 =2j7iJ7iofthe pure Baxter axiom © to derive 

Wf { - k - l W- k) ) = <Pj/ M) J/ M) = \<P{Sf ( - k) ? = l j9(f { - k - 1) ) 2 ^0, 
and again the induction is complete. 



4 The Ring of Methorious Operators 

Let us start by reviewing the general setting for localization in a noncommutative 
unital ring R, following [ 30j § 4.10]. As in the commutative case, the denominator 
set S (the elements that should become invertible) must clearly form a multiplicative 
set S, meaning a submonoid of (R x , •). Clearly we must stipulate 0^5, otherwise the 
localization is the zero ring. If R is a domain, one might want to take S = R x . While 
this is always possible in the commutative setting, one needs an additional condition 
if R is not commutative. 

Indeed, let us strive for a localization S~ l R on the left, meaning all elements have 
the form s~ l r with s £ S and r £ R. Since s , r £ S R this must also be possible for 
the reverse product so that rs~ l — s~ l r for some s £S and f £ R. Multiplying out, we 
get the necessary condition 

SrC\Rs^% for all r £ R and s £ S, (17) 

known as the left Ore condition; the set S is then called left permutable. 

If R is a domain, this condition is actually sufficient for guaranteeing the exis- 
tence of a unique localization S~ l R 3 R, which can be constructed essentially as in 
the commutative case. However, if R has zero divisors, in general one does not get 
an embedding R C S~ l R. In this case, the extension e: R — > S~ l R is a ring homo- 
morphism that is not injective (and of course not surjective). Its kernel contains at 
least those r £ R that yield sr = for some s £ S since this implies e(s)e(r) = and 
hence e(r) = 0, due to e(s) being invertible in S -1 /?. In the classical localization S~ R, 
the kernel should be optimal in the sense that it contains no other elements than these 
necessary ones. 

Definition 6 Let R be an arbitrary ring with SCR. Then e : R — > S~ l R is called a left 

ring of fractions if 

(a) all elements e(s) with s £ S are invertible in S~ l R, 

(b) every element of S~ l R has the form e(s)~ 1 e(r) for some s £ S, r £ R, 

(c) and the kernel of e is given by {r £ R | sr — for some s £ S}. 

The ring homomorphism e is called the extension. 

The missing condition is now easy to establish. Assume s £ S is a right zero 
divisor so that rs — for some r £ R. Then also e(r)e(s) = and hence e(r) = 
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since e(s) is invertible in S l R. But this implies sr = for some s £ S by item (c) of 
Definition|6] Accordingly, one calls a set S with the property 



left reversible. Together with the left Ore condition ( TTTb . this turns out to be sufficient 
for the existence of a left ring of fractions. 

Theorem 2 Let R be an arbitrary ring. Then for any SCR, the left ring of frac- 
tions S~^R exists iffS is multiplicative, left permutable and left reversible. 

Proof See H Thm. 10.6]. 

The setting would become much nicer when the extension e : R — > S _1 R is injec- 
tive so that we can reg ardflCS l R as an embedding. Unfortunately, the localization 
of J£"$[<3, J] that we will work out in the sequel is not of this type. In fact, one can 
easily show that the injectivity of e is equivalent to having only regular elements in S. 
Following 1 9., §5.1], we call an element s £ S regular if it is both left and right regular, 
where left regular means rs = implies r = for all r £ R while right regular means 
sr = implies r = for all r £ R. As ^&[d, j] contains plenty of zero divisors, it is 
hard to achieve a regular denominator set S. 

Of course there are analogous definitions for the right ring of fractions SR~ l , but 
in general the existence of S~ l R does not imply the existence of SR~ l or vice versa, 
and even when both exist they need not be isomorphic. In fact, we shall be dealing 
with a case that is left permutable (Lemma ITTb but not right permutable (Proposi- 
tion 17). 

The localization that we shall construct is based on the monoid <?[<9]$ of regular 
boundary problems over an integro-differential algebra (^,d, J) with umbral char- 
acter set <P. The ring R to be localized is the ring of boundary problems K£\d\&, 
which is a monoid algebra. For such cases, the above setting can be somewhat sim- 
plified. 

A monoid S is called left permutable if it satisfies the Ore condition (T% with 
respect to itself, meaning Ss n SS ^ for all s, S £ S. In the absence of addition, the 
analog of condition (TT~8T > says that for all s, s\ , 52 £ S with s\s — S2S there must exist s £ 
S with ss\ = SS2', in this case we call S left reversible (as a monoid). We call S a left 
Ore monoid if S is both left permutable and left reversible; cf. also Il33l Def. 1.3]. 

The following lemma is a special case of [48, Lem. 6.6] but we include its proof 
since it dispenses with the technical machinery used there. It tells us that we can 
transfer left permutability from S to the monoid algebra KS. 

Lemma 8 IfS is a left permutable monoid, then S is a left permutable subset ofKS. 
Similarly, if S is a left reversible monoid, then S is a left reversible subset ofKS. 

Proof Given s, E S, A, £ K and s £ S, we have to find s, £ S, A,- £ K and s £ S such 



(Vr £ R) (0 £ rS => £ Sr) 



(18) 



that 



s(k\si H h X„s, 



) = (Ajii H hA„s„)s. 



(19) 
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Using the Ore condition in S, we can successively find l n , . . . , l\ s S and r\, . . . ,f n <ES 
such that 

l n S\ = hs 

4-i (/ B *2 ) = hs 
h{h---l n Sn) = r„s 

is fulfilled. Multiplying these equations by Ai l\ ■ ■ ■ l„-\, A2 l\ ■ ■ ■ l n -i, A„ on the 
left yields the system 

s(Alii) = (Xisi)s 
s(X 2 s 2 ) = (X 2 s 2 )s 

5(XnSn) = (KS n )s 

if we set S = l\ ■ ■ ■ l n and Si = h • • • ^j-i'f/ with coefficients A,- = A,-. Summing these 
equations gives the desired Ore condition dl9) . The proof of the second statement 
follows immediately by induction on the number of terms. 

Corollary 1 IfS is a left Ore monoid, the left ring of fractions S~ l (KS) exists. 

Proof Immediate from Lemma[8]and Theorem[2] 

We shall now prove that <?[d]$ is a left Ore monoid. The multiplication of bound- 
ary problems is realized by the semi-direct product ( TTOb of the multiplicative monoid 
<§[d] acting on the additive monoid of boundary spaces. Projecting onto the first fac- 
tor, it is then clear that the monic differential operators of S"[d] must form a left Ore 
monoid. A coefficient algebra with this property shall be called left extensible. 

Example 6 Left permutability requires some common left multiple for given differen- 
tial operators T\,T2& #[<?]. Hence it may seem tempting to restrict ourselves to those 
coefficient algebras (t?,<3) that even have a least common left multiple lclm(7i,72). 
But this setting is not suitable for our purposes: As far as we know, there are only 
two natural examples: The rational functions § — K(x) allow an adaption of the Eu- 
clidean algorithm fl7l p. 23], and of course one may always take & = K. The first 
case yields a differential field, which excludes the existence of an integro-differential 
structure PTT1 p. 518]. The second case restricts us to differential operators with con- 
stant coefficients, which is excessively restrive. 

Fortunately it turns out that we do not need least common left multiples since 
left extensible coefficient algebras are easy to come by, as the next lemma shows. In 
particular, any integro-differential algebra & allows § — K[x] so that <o[d] —A\{K) 
is the Weyl algebra. As another example, consider the standard setting of ExampleQ] 
where for S one may take the larger ring of analytic functions C ffl (R). 

Proposition 6 Any left Noetherian differential domain ((o,d) is left extensible. 
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Proof If § is a left Noetherian domain, it satisfies left permutability by J9] Thm. 5.4] 
and is therefore a left Ore domain. But § [d] is a special case of a skew-polynomial al- 
gebra and therefore inherits the property of being a left Ore domain by [ 9 Prop. 5 .9] . 
But this is not enough since we need the monic differential operators of S"[d] to form 
a left Ore monoid. 

Hence assume T\T = T^T for some monic differential operators T\,T2,T G (o[d]. 
Since <f[<9] is a domain, (Ti — T2)T = implies T\ — T% because T — is not possible. 
Hence the monic differential operators from S"[d] form a left reversible monoid. But 
this monoid is indeed a left Ore monoid: According to [ 10 Lem. 1.5.1], the set of all 
monic polynomials in the skew polynomial ring &[d] is left permutable whenever S 
is a left Noetherian domain. 

Let us remark that the generalization from least common left multiples to com- 
mon left multiples is crucial even for the basic case S — K[x]. For example when T\ = 
d +x and T2 = d 2 +xd +x+ 1, the least common left multiple of T\ and T2 in K{x) [d] 
is 

d 3 + 2 ^d 2 + (x 2 + l+x)d + ^^, 

and so gives T = x<9 3 + (2x 2 — l)d 2 + (x 3 +x 2 +x)d + (x 3 +x 2 — 1) when taken 
in (o[d]. This shows that the least common left multiple of two monic operators 
in (o[d] need not be monic. The conventional computation tools for the Weyl alge- 
bra are therefore not directly applicable for computing monic common left multiples 
in <S\d], but some recent methods [6 53| can be adapted to this purpose [A. Bostan, 
private communication]. 

As explained at the beginning of Section[3] the main tool for ensuring the left Ore 
condition in <§ [d ] is the embedding of singular boundary problems into regular ones . 
This will be achieved in the Regularization Lemma [10] by successively embedding 
the Stieltjes conditions generating a given boundary space. Hence the crucial step is 
to embed a single Stieltjes condition, which we require to be umbral so that we can 
build up a polynomial fundamental system. 

Lemma 9 Let j3 € J^$[<9, j] be an umbral Stieltjes condition over a given integro- 
differential algebra ^ . Then there is a k € N such that (d k+l , [e, . . . , Ed k ~ l , j3]) is a 
regular boundary problem. 

Proof Since /3 is umbral, there is a minimal mononomial x k with jS(x*) 7^ 0. Clearly, 
u = (l,x/l\,--- ,x k /k\) is a fundamental system for d k+ . Using the boundary con- 
ditions 7= (e, . . . , Ed k ~ l we obtain the (« + 1) x (n + 1) evaluation matrix 

/l ••• \ 

Ol'-.i 

7(") - : 

••• 1 

\0 ••• p(x k /k\)J 

where the off-diagonal bottom entries vanish by the assumption on x k . 
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Note that k in Lemma [9] can also be zero, for example if j3 = e. In that case, it 
is not necessary to add any initial conditions since the first-order problem (<9, [e]) is 
already regular. 

Example 7 As a more or less typical case consider /3 = Ei — Eo. Here the minimal 
mononomial is x, which leads back to our standard Example [2] the regular boundary 
problem (d 2 , [eo,Ei — Eo]) = (d 2 , [eq, Ei]). As mentioned earlier, we cannot split off 
the first-order subproblem (d . [e\ — Eo]) since it is singular (this is an example of 
a subproblem that is not a right factor): If it were not, we would know from the 
Division Lemma [3] that the corresponding unique left factor is (<9, [J^]), as shown 
in 1 4 1 . Ex. 28] and also below (Example 8). But of course multiplying out just gives 
the degenerate problem 

(5,[J ]).(5,[e 1 -eo] = (5 2 ,[e 1 -eoD 

and not the desired regularized problem (<3 2 , [eo, Ei]). 

Based on the case of a single boundary condition, we can now embed an arbitrary 
boundary problem into a regular one — provided we work with an umbral character 
set. This is the subject of the following Regularization Lemma. 

Lemma 10 Let <P be an umbral character set for an integro-differential algebra jF. 
Then for an arbitrary boundary problem (T,£t§) £ x K<P there is a regular 
boundary problem (S,£/) G <?[<9],j> that has (T,£$) as a subproblem. 

Proof Let (T,£§) £ K K<P be an arbitrary but fixed boundary problem with T 
a differential operator of order n > and 38 = [j3i , . . . , j3 m ] < of dimension m. 
We write ^ n for the space of initial conditions [e, . . . , Ed" -1 ] and SBk for the partial 
boundary space [jSi , . . . , J3jj with the convention that 3§q = O. We will now prove that 
for each k = 0, . . . , m there is a regular boundary problem (S, £/) that has (T, 3$£) as 
a subproblem. Taking k = m, the theorem follows. 

We use induction on k. The base case follows by setting (S,#f) — (T, J? n )- For 
the induction step, assume is a regular boundary problem that has (T,^_i) 

as a subproblem. We have to construct a regular problem (S,jzf) that has (T,3§k) as 
a subproblem. Letting G be the Green's operator of (§,£/), assume first that faG is 
degenerate. Since /^.G vanishes on J£", we obtain 

Im(G) = J L < [Pk] 1 - and hence [j5 k ] < J, 

so we may set (S,&/) — (§,£/) in that case. Now assume j3^G is nondegenerate and 
hence umbral. Lemma|9]yields a regular problem (T — (d r+l , [e, . . . , E«9 r— 1 , j3^<5]). 
We define the boundary problem 

= (f, ^) ■(§,£/) = (fS, [E§,...,Ed r - l S,p k GS}+s/) 

= (TS, [ES 1 ...,Ed r - x SM+^)i 

where the last equality follows since the conditions /^m = and faGSu = are equiv- 
alent for u 6 sd L . The boundary problem (S, &/) is clearly regular since it is the prod- 
uct of two regular problems, and it has (T, 38^) as a subproblem because 3§k-\ < 
by the induction hypothesis. 
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We are now ready to prove left permutability for the monoid #[<9]$ by merging 
the given factors into a singular problem that is subsequently embedded into a regular 
problem by virtue of Lemma [TOl 

Lemma 11 Let <P be an umbral character set for an integro-differential algebra & 
with left extensible coefficient algebra S. Then #[c?]$ is a left permutable monoid. 

Proof Given {Ti,38i),{T 2 ^%) G <?[<?]<j>, we must find (7i, J^), (f 2 , M 2 ) € g[d]& 
such that (fi, J 1 !) ■{Ti,3§i) = (f 2 , #2) ■ (li,^)- Since S is left extensible, Propo- 
sition [6] yields a common left multiple T and cofactors T\ and ?2 such that T = 
T\ T\ = T 2 T%. Now set SB = 3 \ + S§2- By LemmafTUl there is a regular boundary prob- 
lem (S, sf) that has {T,SS) as a subproblem. But then the boundary problem (S,s</) 
has (7] , ^1 ) and (r 2 , SB2) as regular subproblems, and Lemma[3]yields regular bound- 
ary problems (7i,^i) and (^2,^2) sucn tnat 

(s,*o = Pum) = (72,^2) • (t 2 m 

as claimed. 

Remark 1 Lemma [TTlis also true if one replaces ^[d]^ by the monoid of well-posed 
boundary problems defined before Eq. ( ITTb . This can be seen readily by inspecting 
the above proof (and the proofs of Lemma[3land[T0li. 

Let us give a simple example of a nontrivial Ore quadruple (Ti,S@i), (T2, S§ 2 )> 
(f\ ,S3\), (T 2 ,3@2) that will also serve a good purpose later on. 

Example 8 In the standard setting of ExampleQ] consider the two simplest first-order 
problems on [0,1], namely 

(Zi,#i) = (d,[E<,]) and (T 2 ,m) = (d,[E 1 ]). 

In that case, we have of course T — d, and we apply Lemma [TOl to the boundary 
problem (S,£/) = (d, [eo,Ei]). We end up with 

(a,[ji]).(^[Eo]) = (<9 2 ,[E ,E 1 ]), 

which is already regular. The corresponding cofactors are then 

(T h m) = (d,[!l]) and also (f 2 ,J 2 ) = (d, [J*]) 

since clearly (<3 2 , [j l Q d, e ]) = (<? 2 , [j l Q d, Ei]) = (<? 2 , [e , Ei]). 

As mentioned before, #[t?]$ is left permutable but not right permutable. This 
means there are boundary problems that do not have a common right factor. Actually, 
more is true: Even if we start from two distinct problems with the same differential 
operator, any common right multiple comes from a singular factor. 
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Proposition 7 Let <P be an arbitrary character set for an integro-differential alge- 
bra with coefficient algebra S . Assume (T,Mi), (T,£§2) G have a common 
right multiple 

{T,3§ l ){S&) = {T,3§ 2 ){S^2) (20) 

for some right factors (S,%), (S,^) G £[d] x K4>. Then both (S^i) and {S,^ 2 ) 
are singular whenever @$\ ^ 3%2- 

Proof Assume 8$\ ^ £§2- Projecting onto the boundary spaces, we have BS\S + ^\ = 
@ 2 S + If @\ = O or 3§ 2 = O, we have T = 1 , so in fact 3§i=3§ 2 . Choosing j3i G 
3§\ \3§2 we have j3iS G ^ 2 5+^2 but J3i5 ^ Hence we can write frS = j3 2 S+ y 2 
for some j3 2 £ ^2 and nonzero y 2 G ^2- But then 72 = (J3i - ^2)5 G Ker(5) x n %, 
implies that (S,^) is singular. By symmetry, we see that (S,^) is singular as well. 

Having established that <f[<9]$ is left permutable, the only thing missing for the 
localization is left reversibility, which is very easy in our case. Hence we obtain the 
desired left ring of fractions. 

Theorem 3 Let <t> be an umbral character set for the integro-differential algebra 
(J^\<3, J) with left extensible coefficient algebra S. Then there exists a left ring 
of fractions KS'ld]^, of the ring of boundary problems K$[d]$, with denominator 
set £[d\&. 

Proof By Corollary Q] it suffices to check that is an Ore monoid. Since left 

permutability of has been shown in Lemma fTTl it remains to show that S"[d]<p 

is a left reversible monoid. So assume we have (T\,33\)(T,33) = (T2,3§2)(T,3&) 
for some regular problems [Ti,&i),(T2,$&2),(T,3ff). By left extensibility of § we 
have TT\ — TT2 for some monic T G S\d\ of order n. Let 3$ = [e, • • • , Ed" -1 ] be the 
corresponding space of initial conditions so that (T,0) is a regular problem. Then 
we obtain the regular product problem 

((M)(7i,#i)) (T,&) = ((f , J)(T 2) J)) (T,0), 

where the two parenthesized factors have the same differential operator by the choice 
of T. Since (T, 88) is a regular right factor of this problem, the Division Lemma[3]im- 
plies that (f , 3S)(Ti,3S\) = (f , J , )(7 2 , J 1 ). We conclude that the monic differential 
operators of S"[d] are indeed left reversible. 

The fraction ring KS"[d]^, shall be called the ring of methorious operators (from 
the Greek word [ie-dopLoc meaning "making up a boundary"). Note that it exists in 
particular in the smooth setting: If one chooses J£" = C°°(R), every character is um- 
bral (Proposition H}, and any Noetherian domain can be used for the coefficient al- 
gebra (Proposition^, in particular $ — K[x] or S — C ffl (R). In the sequel, we shall 
confine ourselves to the latter setting. 

As already observed after Theorem|2] we must expect that e : K£[d]<$ — > KS'[d]* ( p 
is not an embedding. For example, we have N = (d, [eo]) — (d, [e\]) G Ker(e) in 
the standard setting of & = C°°(R). For checking this, we use the characterization 
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of Ker(e) given in Definition^)] It suffices to find a regular problem that annihilates N 
from the left. Indeed, we have (d, [J ])N = since we know from Example[8]that 

(d, [Sl})(d, [E ]) = (d\ [E 0lEl ]) = (d, Ul])(d, [ El ]). 

At the moment we do not know Ker(e) in explicit form. However, we have the fol- 
lowing conjecture. 

Conjecture 1 Let <£ be an umbral character set for an integro-differential algebra J? 
with left extensible coefficient algebra S, and let e : K$\d\® — > K^ld]^, be the ex- 
tension into the ring of methorious operators. Then we have A; (7j-,«$?;) 6 Ker(e) 
iff E;A,-G; S where G, is the Green's operator of (T/,^ 1 ,-). 

The previous example is a case in point. The intuitive reason for our conjecture is 
this: Generically, a linear combination of Green's operators has a finite-dimensional 
cokernel (meaning its image is annihilated by just finitely many functionals — in the 
case of a single operator these functionals are the given boundary conditions). Such an 
operator is in some sense "almost invertible". But if a linear combination degenerates 
into an element of (4>), its image becomes one-dimensional, and we cannot expect to 
invert such an "operator". 

5 The Module of Methorious Functions 

Extending an operator ring is much more useful if its elements may still be viewed 
as operating on some — presumably extended — domain of "functions". As explained 
in the Introduction, the operational calculus of Mikusinski avoids this problem by 
merging "operators" and "operands" in the Mikusinski field. In contrast, we shall 
follow the algebraic analysis approach of keeping operators and operands in separate 
structures. 

To this end we construct a suitable module of fractions that extends a given 
integro-differential algebra (^,d, J). In analogy to the Mikusinski approach we 
would have to start from a ring R of integral operators acting on & ', construct its 
localization R* and let it act on the corresponding localization JF*. However, this 
does not work for the following reason: The natural candidate for R would be the 
Zf-algebra of Green's operators for regular boundary problems (since this is the only 
ring over which we have sufficient control). But this monoid/ring is dual to the regu- 
lar boundary problems, as we know from (0. Since is a left Ore monoid but 
not a right Ore monoid, R is a right but not a left Ore ring. Hence we can only create 
a right ring of fractions R* and correspondigly only a right module of fractions that 
should extend & '. But the given integro-differential algebra & only has a natural left 
action of R C ^[d, J]. 

In fact, it is not quite true that one cannot construct a left module of fractions 
from a given left module M over a right Ore ring R. If R* is the right ring of fractions, 
one may of course construct the usual scalar extension M* = R* ®rM. One may then 
refer to M* as a left module of fractions ll48l Prop. 7.2]. But the problem with this 
construction is that we do not know anything about its structure. For example, there 
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is no useful characterization of the kernel of M -+ M* as in the Ore construction of 
Theorem|4]below. 

We must therefore take the differential operators as our starting point, and this 
is why we have constructed the localization KS\d\*^, of the ring of boundary prob- 
lems K£[d]<f,. One might be tempted to take just the simpler ring (o[d] instead of the 
unwieldy monoid ring K$[d\ But this would be too simplistic: In that case one gets 
a two-sided inverse <9~' of the differential operator, no new functions are generated, 
and of course we cannot tackle boundary problems in such a setting. It is essential to 
work with the richer structure K<f>[d]<p, for which we will have to set up a suitable 
left action on # ', 

But first let us briefly review the general setting for the localization of modules. 
As in Theorem^ we go back to the setting of general rings R. In this case we can 
construct the module of fractions in pretty much the same way as the ring of frac- 
tions ED Cor. H.3. 3]. 

Theorem 4 Let M be a left R-module, and let S C R be a multiplicative, right per- 
mutable and right reversible denominator set SCS. Then there exists a left S R- 
module S^^M, and the kernel of the extension jj,: M — > S~ l M consists of those u € M 
for which there exists an s € S with su = 0. 

Now let (J^,<3, J) be an integro-differential algebra and J^"#[<3, J] its ring of op- 
erators. For setting up a suitable left action of ff[<9]0 on J?, we recall the fundamental 
formula OTTl of the Mikusihski calculus 



where s denotes the inverse of the standard convolution operator Z, defined by I f = 
J'q/(<^ ) dt; , and 5q is simply the multiplicative identity generated by the construction 
of the Mikusihski field. Of course one may write 5q = 1 but we keep the explicit no- 
tation for intuitive reasons since we can interpret ( f2TT > as the distributional derivative 
of a function /: 1R — > R that is continuous except for a jump at 0. In that case, 5q is 
the Dirac distribution concentrated at 0. At any rate, the intuitive content of (fJTJ is 
that s is a kind of enhanced differential operator that memorizes the function value 
at that is otherwise lost. Taking this clue, we define the action of ^[dj^ on & by 



where P is the projector onto Ker(T) along S3 L . In the special case of T = d, this 
recovers d2TT) if we think of (d, [eq]) as some kind of algebraic representation of the 
Dirac distribution 8q. 

The definition (l22l i has the consequence that we must extend so that it contains 
the new elements (T,£§), prior to the extension effected by the localization. This is 
somewhat unsatisfying but appears to be inevitable in view of the facts pointed out 
at the opening of the section. Perhaps in the future one can find a more powerful 
localization that generates at once the extension objects needed in the fundamental 
formula (fJTJ. At any rate, we must now also define the action of <£?[<9]$ on the new 
elements, but this is clearly 



sf = f' + f(0)5o, 



(21) 



(T,&)-f = Tf + Pf(T,&) 



(22) 



{f,&)-f{T,38)=f{TT,$T + 38), 



(23) 
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where the boundary problem on the right side is the usual product ( TTOb in the problem 
monoid <§ [d]<p. This will ensure that the action is well-defined. 

The intuitive meaning of the ideal element f{T,SB) is that it records various 
integration constants that were lost under differentation. When the element is created 
by (|22| |. it encodes all the boundary values j3 (/) for j3 € SB, and the function / is 
confined to the kernel of T. Further action by a boundary problem (T,SB) according 
to d23l ) leads to a new differential operator TT while leaving / intact. Hence we see 
that the function inside an ideal element is always in the kernel of the associated 
differential operator. We may therefore restrict the ideal elements / {T,SB) to those 
with Tf = 0. After the action of d23l ). the new boundary space is SBT + SB, so the 
old boundary conditions /3 £ SB are retained (the ones for which we have boundary 
values). The new boundary conditions j5T for j5 e SB give zero on / since actually 
Tf = 0. In some sense, we have added redundant boundary data. 

We must therefore regulate redundancy in the ideal elements. This can be achieved 
by declaring that / (T, SB) should be the same as Gf {T, 3B)(T,SB) for any (f , =#) 6 
<o°[<3]$ with Green's operator G. This can be understood as follows. The boundary 
conditions in Gf (T, SB) (f , SB) are j3 £ J 1 and j3 T with J3 £ SB. The former yield zero 
on Gf since G is the Green's operator of (T \SB), the latter give back j3 f(Gf) = j3(/) 
just as for the ideal element / (T, SB). This leads to the following definition. 

Definition 7 Let ,f be the subspace of & ®k K$[d\& generated by the ideal ele- 
ments f®{T,SB)=f(T,SB) with Tf = 0. Furthermore, let ,y Q be the subspace of J 
generated by the elements 

/ (r, SB) - Gf (Tf,SBf + SB). 
Then we define the module of methorious functions to be = & ®J?/J?o- 
Let us now check that is indeed a module. 

Proposition 8 Let (J?, <9, J) be an integro-differential algebra with character set <t>. 
The definitions (122b and d231 > induce a monoid action of $[d]$> on such that it 
becomes a KS'[d]<p-module. 

Proof Any monoid E acting on a Zf-algebra A extends to an action of the monoid 
ring K[E] that makes A into a K[E] -module. Hence it suffices to verify the statement 
about the monoid action. 

First of all we must check that (l23l l does not depend on the representative / (7*, SB) € 
J?, meaning it maps the subspace J^o into itself. This follows immediately from the 
associativity of the multiplication ( TTOb in 



Clearly we have (1 , 0) • / = / + (1 , 0) = / and (1 , 0) • (T, SB) f = (T, SB) f, so the 
unit element is respected. Associativity is immediate in the case of d23T >. so it remains 
to check that {T, SB) ■ ((? ,SB) ■ f) and ((T, SB)(T,SB)) ■ f are equal. The former is 



where P is the projector onto Ker(T) along SB and P is correspondingly the projector 
onto Ker(T') along If G and G are the Green's operators of (T,SB) and (T,3&), 
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respectively, one can easily check that the projector associated with the composite 
problem 

(T,M)(f , S3) = {Tf,3$f + &) 

is P + GPT, meaning it projects onto Ker(7T) along (SSt + £$)^. Hence we obtain 
for the other side of the prospective equality 

(TT,3§f + M)-f = TTf + {Pf + GPtf) (TT,&T + 33) 
= TTf + GPTf(TT,3§T + 3§)+Pf(TT,3§T + 3§), 

where the first and the last term is identical to the corresponding terms on the left 
hand side while the middle terms are equal since their difference is in J?q. 

Since is a ^^[dj^-module, we obtain the KS'ld]^ -module J^, of metho- 
rious hype rf unctions by localization via Theorem [4] We cannot expect the exten- 
sion n : -J^<p — > J^J to be injective since its kernel 

jf = {<pe | 3(T,@) e£[d]&: (r,^) • <p = 0} (24) 

contains elements like (<9,[ei]) — (<9,[eo]), which is annihilated upon multiplying 
with (<5,[Jq]) from the left (see the example before Conjecture [TJ. But fortunately 
no elements of are lost. 

Proposition 9 Let <9, J) be an integro-differential algebra with character set <P. 
Then we have an embedding C 

Proof Assume {T,S§)-f = for some / G & and (T,3§) G g[d]&. Then Tf + 
Pf (T, SS) = implies that Tf = and Pf = 0. But the former means that / G Ker(r) 
and the latter that / G S3 L . Since (T,3§) is a regular problem, we have a direct 
sum Ker(r) + 3g x = & and so / = 0. 

In the module of methorious hyperfunctions, we can finally justify our earlier 
notation (T, £8) ~ 1 for the Green 's operator of a regular boundary problem (T, 33) G 
£"[d}&. To avoid confusion with (T,3§)~ 1 G KS\d\^ we refrain from the reciprocal 
notation for Green's operators in the scope of the following theorem. But the result of 
the theorem is of course that it anyway does not matter how we interpret (T,3§)~ l f 
since it amounts to the same. 

Proposition 10 We have (T \3§)~ l ■ f = Gf for all f G Moreover, if T f = 
then (J,2$Y X -f(T,M)=f. 

Proof We have (T,38) ■ Gf = TGf + PGf {T,3§) = f since TG = 1 and PG = 0. 
Multiplying by (T, 3§)~ l G /f<?[<9]^ yields the first result claimed. Now assume Tf = 
0. We obtain (T,3§) ■ f = + Pf (T,&) = f (T, S3) since Pf = / for / G Ker(T). 
Again we multiply by (T,3§)~ 1 to gain the result. 

We conclude with an example that hints at possible applications of our noncom- 
mutative Mikusihski calculus. The classical Mikusiriski calculus has only one funda- 
mental formula since boundary values (or rather: initial values) are only processed 
at 0. In contrast, there are plenty of fundamental formulae in the noncommutative 
Mikusinski calculus. 
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Example 9 Let us write and 8^ as abbreviations for the problems (<9,[e,e]) G 
KS[d]% and (d, [e § ]) G J^, respectively. Then we have ^/ = /' + /(§) 5§ by 
the definition of the action d22l) . So we have algebraic representations for all the 
Dirac distributions. But there are other methorious functions that do not have any 
distributional counterpart. For example, let us consider df = (<3,[J ]) G KS[d\^ 
and e = (d, [J ]) G This yields the fundamental formula 

d F f = f'+(j^m)d^je, 

so £ is a kind of "smeared out" distribution that keeps the mean value of a given 
function /. 

Example 10 Finally let us see how one solves inhomogeneous boundary problems in 
the noncommutative Mikusinski calculus. We use the standard setting of = C°°(]R) 
and consider the problem 



u(0) = a,u(b) = b ' 

for given boundary values a,b£ R. We know how to compute the Green's operator G 
of the homogeneous problem (d 2 , [eo,Ei]) as mentioned at the end of Section [2] 
In this case the projector P onto Ker(<9 2 ) = along [eojEi]^ is given by Pu = 
k(0) (1 — x) + u(l)x. Hence we have 

(d 2 , [Eo.Ei]) U = ll" + (ll(0) (1 -X) + Il(l)jf) (d\ [Eo,Ei]) 

= f+(a(l-x)+b x yd 2 ,[E ,E l }), 
u = Gf + a (1 — x)+i>x 

by PropositionfTOl 



6 Conclusion 

As already indicated at various points above, our construction has several loose ends. 
To begin with, it would be preferrable to localize in a subring (or even all) of [d , J] 
rather than the monoid ring K$[d]®. This would have the advantage that one has a 
natural action on the underlying integro-differential algebra (JF,d,j), and the some- 
what artificial action on the module of methorious functions would be unnec- 
essary. Regarding the latter, we have already remarked in Section [5] that the current 
two-stage process of creating the localization J£"J is unsatisfactory: In Mikusiriski's 
setup, all the "ideal elements" like s and 5o are an immediate result of the localization 
while we have to supply them offhand, prior to localization. We would like to find a 
better formulation in the future that will avoid this kind of inadequacy. 

One possible path towards such an improved localization is suggested by the 
prominent appearance of singular problems in the proof of Lemma QT| It may be 
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worthwhile to expand the ring KS[d\& to include (all or some) singular boundary 
problems. Of course this will also increase the kernel of the corresponding exten- 
sion (we cannot expect to invert singular problems), but perhaps the resulting ring 
of fractions is more natural. In particular, it may be possible to interpret the sum of 
two (singular) boundary problems in some useful way, in contrast to the formal sums 
of K£[d\&. The results of 02711251 will be useful for the work in this direction. 

Staying with the current construction, there are some obvious open questions: 
Fist of all, the kernel of the extension £ into the methorious operators should be deter- 
mined, possibly by proving CojectureQ] Likewise, the kernel [24] of the extension /i 
into the methorious hyperfunctions is to be computed. 

Finally, we would like to draw attention to the interesting link between integro- 
differential algebras and the umbral calculus (Sectional, which deserves to be studied 
in more depth. There is also an intriguing relation [ 19] between the umbral calculus 
and the Rota-Baxter algebras (both of these being favourite topic of G.-C. Rota), 
which might benefit from the new perspective afforded by boundary problems in 
integro-differential algebras (since the latter are special cases of differential Rota- 
Baxter algebras and hence of plain Rota-Baxter algebras). As a more mundane goal, it 
would also be important to find a better characterization of umbral character sets that 
strengthens the separativity and completeness conditions given before Example]!] 

Acknowledgements First of all we would like to acknowledge Georg Regensburger, who has worked with 
the first author on an earlier (much more specialized and mathematically flawed) version of localization. 
We thank him also for his insight about monic common left multiples and the corresponding reference 
(last paragraph of Proposition [6). We express our gratitude also to A. Bostan for his pertinent remarks on 
left common multiples in the Weyl algebra (see the remark after the proof of Proposition[6). Furthermore, 
we are grateful for illuminating discussions with I. Dimovski and M. Spiridonova. 



References 

1. Agarwal, R.: Boundary value problems for higher order differential equations. World Scientific Pub- 
lishing Co., Teaneck, NJ (1986) 

2. Agarwal, R., O'Regan, D.: An introduction to ordinary differential equations. Springer (2008) 

3. Albrecher, H., Constantinescu, C, Pirsic, G., Regensburger, G., Rosenkranz, M.: An algebraic ap- 
proach to the analysis of gerber-shiu functions. Insurance: Mathematics and Economics Special Issue 
on Gerber-Shiu Functions, accpeted (2009) 

4. Baxter, G.: An analytic problem whose solution follows from a simple algebraic identity. Pacific J. 
Math. 10, 731-742 (1960) 

5. Berg, L.: Einfiihrung in die Operatorenrechung. Deutscher Verlag der Wissenschaften, Berlin (1962) 

6. Bostan, A., Chyzak, E, Li, Z., Salvy, B.: Fast computation of common left multiples of linear ordinary 
differential operators. ACM Commun. Com put. Algebra 45(1/2), 111-112 (201 1). DOI 10.1145/ 
2016567.2016581. URL http://doi.acm.org/10. 1145/2016567.2016581 

7. Bourbaki, N.: Topological vector spaces. Chapters 1-5. Elements of Mathematics (Berlin). Springer, 
Berlin (1987) 

8. Chyzak, F.: Holonomic systems and automatic proofs of identities. Tech. Rep. 2371, Institut National 
de la Recherche en Informatique et Automatique (1994) 

9. Cohn, P.M.: Introduction to Ring Theory. Springer (2000) 

10. Cohn, P.M.: Free ideal rings and localization in general rings, New Mathematical Monographs, 
vol. 3. Ca mbridge Universi ty Press, Ca mbridge (2006). DOI 10.1017/CBO9780511542794. URL 
|http : //dx ■ dol ■ org/10 ■ 1017 /CB0978051 1542794] 

11. DiBucchianico, A.: An introduction to umbral calculus (1998). Lectue notes, available at 
www . win . tue . nl/~ adibucch/eidma . ps 



A Noncommutative Mikusinski Calculus 



31 



12. Dimovski, I.: Convolutional calculus, Mathematics and its Applications (East European Series), 
vol. 43. Kluwer Academic Publishers Group, Dordrecht (1990) 

13. Dimovski, I.: Nonlocal operational calculi. Trudy Mat. Inst. Steklov. Izbran. Voprosy Mat. Fiz. i 
Anal. 203, 58-73 (1994) 

14. Dimovski, I.: Operational calculi for boundary value problems (2012). Talk at the Algebraic Aspects 
of Differential and Integral Operators Session (AADIOS'12), 18th Conference on Applications of 
Computer Algebra (ACA), Sofia, 25-28 June 

15. Duffy, D.G.: Green's functions with applications. Studies in Advanced Mathematics. Chapman & 
Hall, Boca Raton, FL (2001) 

16. Engl, H.W.: Integralgleichungen. Springer Lehrbuch Mathematik. Springer- Verlag, Vienna (1997) 

17. Flowe, R.P., Harris, G.A.: A note on generalized Vandermonde determinants. Siam J. Matrix anal, 
appl 14, 1146-1151 (1993) 

18. Grigoriev, D.Y.: Complexity of factoring and calculating the GCD of linear ordinary differential op- 
erators. J. Symbolic Comput. 10(1), 7-37 (1990) 

19. Guo, L.: Baxter algebras and the umbral calculus. Adv. in A ppl. Math. 27(2- 3), 405-426 (2001). 
DOI 10.1006/aama.2001.0742. URL |http : //dx . doi . o rg/10 . 1006/aama . 20 01 .0742 1 Special 
issue in honor of Dominique Foata's 65th birthday (Philadelphia, PA, 2000) 

20. Guo, L.: Baxter algebras and differential algebras. In: Differential algebra and related topics (Newark, 
NJ, 2000), pp. 281-305. World Sci. Publ., River Edge, NJ (2002) 

21. Heaviside, O.: On operators in physical mathematics, part i. In: Proceedings of the Royal Society, 
vol. 52, pp. 504-529. London (1893) 

22. Heaviside, O.: On operators in physical mathematics, part ii. In: Proceedings of the Royal Society, 
vol. 54, pp. 105-143. London (1894) 

23. Kamke, E.: Differentialgleichungen. Lbsungsmethoden und Losungen. Teil I: Gewohnliche Differ- 
entialgleichungen, Mathematik und ihre Anwendungen in Physik und Technik A, vol. 18, eighth edn. 
Akademische Verlagsgesellschaft, Leipzig (1967) 

24. Kolchin, E.: Differential algebra and algebraic groups, Pure and Applied Mathematics, vol. 54. Aca- 
demic Press, New York- London (1973) 

25. Korporal, A.: Symbolic methods for generalized green's operators and boundary problems. Ph.D. 
thesis, Johannes Kepler University, Linz, Austria (2012). In progress 

26. Korporal, A., Regensburger, G, Rosenkranz, M.: A Maple package for integro-differential operators 
and boundary problems (2010). Also presented as a poster at ISSAC ' 10. 

27. Korporal, A., Regensburger, G, Rosenkranz, M.: Regular and singular boundary problems in maple. 
In: Proceedings of the 13th International Workshop on Computer Algebra in Scientific Computing, 
CASC'201 1 (Kassel, Germany, September 5-9, 201 1), Lecture Notes in Computer Science, vol. 6885. 
Springer (2011) 

28. Korporal, A., Regensburger, G, Rosenkranz, M.: Symbolic computation for ordinary boundary prob- 
lems in maple. In: Proceedings of the 37th International Symposium on Symbolic and Algebraic 
Computation (ISSAC'09) (2012). Software presentation 

29. Krattenthaler, C: Advanced determinant calculus. Sem. Lothar. Combin. 42, Art. B42q, 67 pp. (elec- 
tronic) (1999). The Andrews Festschrift (Maratea, 1998) 

30. Lam, T.Y.: Lectures on modules and rings, Graduate Texts in Mathematics, vol. 189. Springer- Verlag, 
New York (1999) 

31. Mikusinski, J.: Operational calculus, International Series of Monographs on Pure and Applied Math- 
ematics, vol. 8. Pergamon Press, New York (1959) 

32. Nachbin, L.: Sur les algebres denses de fonctions differentiables sur une variete. C. R. Acad. Sci. 228, 
1549-1551 (1949) 

33. Picavet, G: Localization with respect to endomorphisms. Semigroup Forum 67(1), 76-96 (2003). 
DOI 10.1007/s00233-002-0008-2. URL |http : //dx . doi . org/10 . 1007/s00233-002-0008-2l 

34. van der Put, M., Singer, M.F: Galois theory of linear differential equations, Grundlehren der Mathe- 
matischen Wissenschaften, vol. 328. Springer, Berlin (2003) 

35. Regensburger, G, Rosenkranz, M.: An algebraic foundation for factoring linear boundary prob- 
lems. Ann. Mat. Pura Appl. (4) 188(1), 123-151 (2009). DOI 10.1007/sl0231-008-0068-3. URL 
|http : //dx . doi . org/10 ■ i007/sl0231-008-0068-3| DOI:10.1007/sl0231-008-0068-3 

36. Regensburger, G, Rosenkranz, M., Middeke, J.: A skew polynomial approach to integro-differential 
operators. In: ISSAC 2009 — Proceedings of the 2009 International Symposium on Symbolic and 
Algebraic Computation, pp. 287-294. ACM, New York (2009) 

37. Ritt, J.F: Differential algebra. Dover Publications Inc., New York (1966) 



32 



Markus Rosenkranz, Anja Korporal 



38. Rosenkranz, M.: Lagrange Inversion. Master's thesis, Research Institute for Sym- 
bolic Computation (RISC), Johannes Kepler University, Linz (1997). Available at 
http : //www . rise . jku. at/publications/download/risc_2383/Diploma.pdf 

39. Rosenkranz, M.: A new symbolic method for solving linear two-point boundary value problems on 
the level of operators. J. Symbolic Comput. 39(2), 171-199 (2005) 

40. Rosenkranz, M., Regensburger, G.: Integro-differential polynomials and operators. In: D. Jeffrey (ed.) 
ISSAC'08: Proceedings of the 2008 International Symposium on Symbolic and Algebraic Computa- 
tion. ACM Press (2008) 

41. Rosenkranz, M., Regensburger, G.: Solving and factoring boundary problems for linear ordinary dif- 
ferential equations in differential algebras. Journal of Symbolic Computation 43(8), 515-544 (2008). 
DOI 10.1016/j.jsc.2007.11.007 

42. Rosenkranz, M., Regensburger, G., Tec, L., Buchberger, B.: A symbolic framework for operations 
on linear boundary problems. In: V.P. Gerdt, E.W. Mayr, E.H. Vorozhtsov (eds.) Computer Algebra 
in Scientific Computing. Proceedings of the 11th International Workshop (CASC 2009), LNCS, vol. 
5743, pp. 269-283. Springer, Berlin (2009) 

43. Rosenkranz, M., Regensburger, G., Tec, L., Buchberger, B.: Symbolic analysis of boundary problems: 
From rewriting to parametrized Grobner bases. In: U. Langer, P. Paule (eds.) Numerical and Symbolic 
Scientific Computing: Progress and Prospects, pp. 273-331. Springer (2012) 

44. Rota, G.C.: Baxter algebras and combinatorial identities (I, II). Bull. Amer. Math. Soc. 75, 325-334 
(1969) 

45. Salvy, B., Zimmerman, P.: Gfun: a maple package for the manipulation of generating and holonomic 
functions in one variable. ACM Trans. Math. Softw. 20(2), 163-177 (1994). DOI http://doi.acm.org/ 
10.1145/178365.178368 

46. Schwarz, R: A factorization algorithm for linear ordinary differential equations. In: ISSAC '89: Pro- 
ceedings of the ACM-SIGSAM 1989 international symposium on Symbolic and algebraic computa- 
tion, pp. 17-25. ACM Press, New York, NY, USA (1989). DOI http://doi.acm.org/10.1145/74540. 
74544 

47. Schwarz, F.: Algorithmic Lie theory for solving ordinary differential equations, Pure and Applied 
Mathematics (Boca Raton), vol. 291. Chapman & Hall/CRC, Boca Raton, FL (2008) 

48. Skoda, Z.: Noncommutative localization in noncommutative geometry. In: Non-commutative lo- 
calization in algebra and topology, London Math. Soc. Lecture Note Sen, vol. 330, pp. 220- 
313. Ca mbridge Univ. Press, Cambridge (2006). DOI 10. 1017/CBO9780511526381.015. URL 
|http : //dx . doi . org/ 10 . 10i7 /CB0978051 1526 381 ■ 015 | 

49. Spiridonova, M.: Operational methods in the environment of a computer algebra system (2010). Talk 
at the Algebraic Aspects of Differential and Integral Operators Session (AADIOS' 10), 16th Confer- 
ence on Applications of Computer Algebra (ACA), Vlora, 24-27 June 

50. Stakgold, I.: Green's functions and boundary value problems. John Wiley & Sons, New York (1979) 

51. Stenstrom, B.: Rings of quotients. Springer- Verlag, New York (1975). Die Grundlehren der Mathe- 
matischen Wissenschaften, Band 217, An introduction to methods of ring theory 

52. Tec, L.: A symbolic framework for general polynomial domains in theorema: Applications to bound- 
ary problems. Ph.D. thesis, Research Institute for Symbolic Computation, Johannes Kepler University, 
A-4040 Linz, Austria (201 1) 

53. Tsai, H.: Algorithms for algebraic analysis. Ph.D. thesis, University of California at Berkeley (2000) 

54. Tsarev, S.P.: An algorithm for complete enumeration of all factorizations of a linear ordinary dif- 
ferential operator. In: ISSAC '96: Proceedings of the 1996 international symposium on Sym- 
bolic and algebraic computation, pp. 226-231. ACM Press, New York, NY, USA (1996). DOI 
http://doi.acm.org/10.1145/236869.237079 



